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Abstract: In this paper we propose a distributed dynamic controller for sharing frequency control 
reserves of asynchronous AC systems connected through a multi-terminal HVDC (MTDC) grid. 
We derive sufficient stability conditions, which guarantee that the frequencies of the AC systems 
converge to the nominal frequency. Simultaneously, the global quadratic cost of power generation 
is minimized, resulting in an optimal distribution of generation control reserves. The proposed 
controller also regulates the voltages of the MTDC grid, asymptotically minimizing a quadratic 
cost function of the deviations from the nominal voltages. The proposed controller is tested on 
a high-order dynamic model of a power system consisting of asynchronous AC grids, modelled 
as IEEE 14 bus networks, connected through a six-terminal HVDC grid. The performance of 
the controller is successfully evaluated through simulation. 


1. INTRODUCTION 

Power transmission over long distances with low losses is 
one of the main challenges in today’s power transmission 
systems. As the share of renewables rises, so does the need 
to balance generation and consumption mismatches, often 
over large geographical areas. Due to the high resistive 
losses in AC cables, high-voltage direct current (HVDC) 
power transmission is a commonly used technology for 
long-distance power transmission. The higher investment 
cost of an HVDC transmission system compared to an AC 
transmission system is compensated by the lower resistive 
losses for sufficiently long distances [Melhem 2013]. The 
break-even point, i.e., when the total construction and 
operation costs of overhead HVDC and AC lines equal, 
is typically 500-800 km [Padiyar 1990]. However, for 
cables, the break-even point is typically less than 50 
km [Van Hertem et al. 2010]. Increased use of HVDC 
technologies for electrical power transmission suggests 
that future HVDC transmission systems are likely to 
consist of multiple terminals connected by several HVDC 
transmission lines [Haileselassie and Uhlen 2013]. Such 
systems are referred to as multi-terminal HVDC (MTDC) 
systems in the literature. Many existing AC grids are 
connected through legacy HVDC links, which are typically 
used for bulk power transfer between AC areas, rather 
than balancing generation mismatches. The fast operation 
of the DC converters however also enables frequency 
regulation of one of the connected AC grids through the 
HVDC link. One example is the frequency regulation of 
the island of Gotland in Sweden, which is connected to the 
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Nordic grid through an HVDC cable [Axelsson et al. 2001]. 
However, since the Nordic grid has orders of magnitudes 
higher inertia than the AC grid of Gotland, the influence 
of the frequency regulation on the Nordic grid is negligible. 

By connecting multiple AC grids by an MTDC system, 
the frequency regulation reserves in each AC grid can 
be shared with connected AC grids, which reduces the 
need for frequency regulation reserves in the individual 
AC systems [Li et al. 2008]. In recent years, this idea 
has gained increasing interest in the literature. Dai et al. 
[2011] and Silva et al. [2012] employ several controllers 
with decentralized structure to share frequency control 
reserves. In [Silva et al. 2012] no stability analysis of the 
closed-loop system is performed, whereas [Dai et al. 2011] 
guarantees stability provided that the connected AC areas 
have identical parameters and the voltage dynamics of 
the HVDC system are neglected. In [Taylor and Scardovi 
2014], an optimal decentralized controller for AC systems 
connected by HVDC systems is derived. In contrast to 
the aforementioned references, Andreasson et al. [2015a] 
consider the dynamics of connected AC systems as well 
as the dynamics of the MTDC system. By connecting 
the AC areas with a communication network, the per¬ 
formance of the controller can be further improved. Dai 
et al. [2010] consider a distributed controller relying on 
a communication network to share frequency control re¬ 
serves of asynchronous AC transmission systems connected 
through an MTDC system. However, the controller re¬ 
quires a slack bus to control the DC voltage, and is thus 
only able to share the generation reserves of the non-slack 
AC areas. Another distributed controller is proposed by 
Dai and Damm [2011]. Stability is guaranteed, and the 
need for a slack bus is eliminated. The voltage dynamics 
of the MTDC system are however neglected. Moreover the 
implementation of the controller requires every controller 
to access measurements of the DC voltages of all other 



MTDC terminals. [Andreasson et al. 2015b] propose a 
distributed secondary generation controller. In contrast to 
the aforementioned references, the MTDC dynamics are 
explicitly modeled, and the DC voltages are controlled in 
addition to the frequencies. 

Despite the coordination through a communication net¬ 
work, the distributed controller by Andreasson et al. 
[2015b] fails in eliminating the static errors of the frequen¬ 
cies. In this paper we address this issue by introducing 
a distributed secondary controller also for the power fed 
into the MTDC grid from the AC grids. This allows us 
to eliminate any static errors in the frequency deviations. 
Furthermore, quadratic cost functions of the power gener¬ 
ation and the voltage deviations are minimized. 

The remainder of this paper is organized as follows. In 
Section 2, the system model and control objectives are 
defined. In Section 3, a distributed secondary frequency 
controller for sharing frequency control and restoration 
reserves is presented, and is shown to satisfy the control 
objectives. In Section 4, simulations of the controller on a 
six-terminal MTDC test system are provided, showing the 
effectiveness of the proposed controller. The paper ends 
with concluding remarks in Section 5. 

2. MODEL AND PROBLEM SETUP 

2.1 Notation 

Let be a static, undirected graph. Denote by V = 
{1,..., n} the vertex set of Q, and by £ = {1,..., m} the 
edge set of Q. Let A/] be the set of neighboring vertices to 
i G V. Denote by B the vertex-edge incidence matrix of Q, 
and let £vv = BWB'^ be the weighted Laplacian matrix of 
Q, with edge-weights given by the elements of the diagonal 
matrix W. We denote the space of real-valued n x m- 
valued matrices by Let C“ denote the open left half 

complex plane, and C“ its closure. We denote by c„xm a 
vector or matrix of dimension n x m, whose elements are 
all equal to c. Eor simplifying notation, we write c„ for 
Cnxi- For a symmetric matrix A, A > 0 (A > 0) is used 
to denote that A is positive (semi) definite. /„ denotes the 
identity matrix of dimension n. For simplicity, we will often 
drop the notion of time dependence of variables, i.e., x{t) 
will be denoted x. Let j]-]]^ denote the maximal absolute 
value of the elements of a vector. 

2.2 Model 

We will give here a unified model for an MTDC sys¬ 
tem interconnected with several mutually asynchronous 
AC systems. We consider an MTDC transmission system 
consisting of n converters, denoted i = l,...,n, each 
connected to an AC system, i.e., there are no pure DC 
nodes of the MTDC grid. The converters are assumed to be 
connected by an MTDC transmission grid. The dynamics 
of converter i is assumed to be given by 

= ( 1 ) 

where Vi is the voltage of converter i, > 0 is its capac¬ 
itance, and is the injected current from an AC grid 
connected to the DC converter. The constant Rij denotes 


the resistance of the HVDC transmission line connecting 
the converters i and j. The graph corresponding to the 
HVDC line connections is assumed to be connected. The 
AC system is assumed to consist of a single generator 
which is connected to the corresponding DC converter, 
representing an aggregate model of the AC grid. The dy¬ 
namics of the AC system are given by the swing equation 
[Machowski et al. 2008]: 

m,uj, = pr+pr-pr, ( 2 ) 

where rrii > 0 is its moment of inertia. The constant 
is the generated power, P™ is the power load and P^-^ is 
the power injected to the DC system through converter i, 
respectively. The powers are all assumed to be deviations 
from a nominal operation point. The control objective can 
now be stated as follows. 

Objective 1. The frequency deviations are asymptotically 
equal to zero, i.e., 

lim uji{t) — = 0 i = 1,..., n, (3) 

where is the nominal frequency. The total quadratic 
cost of the power generation is minimized asymptotically, 
i.e., limt^oo Pf™ = Pf®*"*, Vi = 1,..., n, where 

n 

[Pr *, •. •, Pn = argmin - ^ /f P^ (4) 

i-j 5 n 

subject to Pf®” -I- P™ — P™^ = 0, Vi = l,...,n and 
^"^1 P™-^ = 0, i.e., power balance both in the AC grids 
and in the MTDC grid, in the absence of power losses. 
The positive constants ff represent the local cost of 
generating power. Finally, the voltages are such that the 
a quadratic cost function of the voltage deviations is 
minimized asymptotically, i.e., limt_>oo h) = L)*,Vi = 
1,..., n, where 

n 

[u;,..., IP] = argmin - ^ {V - (5) 

subject to (3)-(4), and where the fY is a positive constant 
reflecting the local cost of voltage deviations, and U/®^ is 
the nominal voltage of converter i. 

Remark 1. The minimization of (4) is equivalent to power 
sharing, where the generated power of AC area i is 
asymptotically inverse proportional to the cost ff. The 
cost ff can be chosen to reflect the available generation 
capacity of area i. 

Remark 2. It is in general not possible that limt_>oo = 
U/®^ Vi = 1,..., n, since this does not allow for the currents 
between the HVDC converters to change by (1). Note that 
the optimal solution to (4) fixes the relative voltages, leav¬ 
ing only the ground voltage as a decision variable of (5). 
Note also that the reference voltages U/®^, i = 1,..., n, are 
generally not uniform, as is the reference frequency w®®^. 
Remark 3. Note that Objective 1 does not include con¬ 
straints of, e.g., generation and line capacities. This re¬ 
quires that the perturbations from the operating point are 
sufficiently small, to guarantee that these constraints are 
not violated. In case of a large disturbance, e.g., a fault, 
the constraints might be violated. 



3. COORDINATED SECONDARY EREQUENCY 
CONTROL 


3.1 Controller structure 

In this section we propose a distributed secondary fre¬ 
quency controller. In addition to the generation controller 
proposed by Andreasson et al. [2015b], we also propose 
a secondary controller for the voltage injections into the 
HVDC grid. The distributed generation controller of the 
AC systems, which was first given by Andreasson et al. 
[2015b], is given by 

Pf" = % 

i)i = c'ljiri, ( 6 ) 

where iLf °°p and A:f™°P’ ^ are positive controller parame¬ 
ters. Moreover c’f = > 0, i.e., the communication graph 

is undirected. The above controller can be interpreted as 
a distributed Pl-controller, with a distributed consensus 
filter acting on the integral states iji. The proposed con¬ 
verter controllers governing the incremental power injec¬ 
tions from the AC systems into the HVDC grid are given 
by 

Pr = + KY{Vr^ - Vi) 

+ H - (i>j) 

j&Vi 

K‘^ 

(l>Z = ( 7 ) 

where and K‘f are positive controller parameters, 
7 > 0 and > 0. The converter controller (7) can 

be interpreted as an emulation of an AC network between 
the isolated AC areas. The auxiliary controller variables 
(j)i can be thought of as the phase angles of AC area z, 
governing the power transfer between the areas, if these 
were connected by AC lines rather than an MTDC grid. 
In contrast to an AC system, however, the power is fed 
into the HVDC grid instead of being directly transferred 
to the connected AC areas. 

The HVDC converter is assumed to be perfect and instan¬ 
taneous, i.e., injected power on the AC side is immediately 
and losslessly converted to DC power. Furthermore, the 
dynamics of the converter are ignored, implying that the 
converter is assumed to track the output of controller 
(7) instantaneously. This assumption is reasonable due 
to the dynamics of the converter typically being orders 
of magnitudes faster than the primary frequency control 
dynamics of the AC system [Kundur 1994]. The relation 
between the injected HVDC current and the injected AC 
power is thus given by 

Vzir = pr. (8) 

By assuming Vi = I/"°™ Vi = 1,... ,n, where I/”°™ is a 
global nominal voltage, we obtain 

vnom/'nj = (9) 

3.2 Stability analysis 

We are now ready to analyze the stability of the closed- 
loop system. Dehne the stacked state vectors as ui = uj — 


and V = V — where w = [wi,... ,a;„]’^, V = 
[V,...,Vnf, w®f = [yr',---,KlT, V = 
and (j) = [</>!,..., (/)„]. Combining the voltage dynamics 
(1), the frequency dynamics (2) and the generation control 
(6), the converter controller (7) and the power-current 
relationship (9), we obtain the closed-loop dynamics 

iij=M(^- + K‘^)oj + K^V 

z 1 / \ 1 

V=- - EK‘^C -E{£r + -- V + -- E£^(j) 

^nom \ ^ p^nom / ^nom 

f] =A:droop.I^ _ 

(10) 

where M = diag(mi“^,... ,mn~^) is a matrix of inverse 
generator inertia, E = diag(C']” ,..., C~^) is a matrix of 
electrical elastances, Er is the weighted Laplacian matrix 
of the MTDC grid with edge-weights IfRij, £r^ and £(f, 
are the weighted Laplacian matrices of the communication 
graphs with edge-weights and c^, respectively, and 

pm _ [p™^... ^ p™]'r_ yye define the diagonal matrices 
of the controller gains as K'^ = diag(Ar“,..., Al])'), etc. 

Let y = [w^, V^]^ define the output of (10). Clearly the 
linear combination 1^(/) is unobservable and marginally 
stable with respect to the dynamics (10), as it lies in the 
nullspace of £^. In order to facilitate the stability analysis, 
we will perform a state-transformation to this unobserv¬ 
able mode. Consider the following state-transformation: 



where S is an n x (n — 1) matrix such that S'I is 

orthonormal. By applying the state-transformation (^11) 
to (10), we obtain dynamics where it can be shown that 
the state </>) is unobservable with respect to the defined 
output. Hence, omitting (])] does not affect the output 
dynamics. Thus, we define 0" = [02) ■ • ■) 0n]) and obtain 
the dynamics 

Cj =m( - (a:'^™°p -p K^)Cj + K^V 

-a:^(a:“)-ia:'^™°P’ - E^Scj)" + p™) 

Z 1 ( K'^ \ .. 1 

V =- - EK^C -E £r+ - - V + - - E£^S(/>'' 

p^nom I ^ ^nom / p^nom ^ ^ 

z) - £^77 

0"=S'^(iVP’)-iA:“w-70". (12) 

Provided that the system matrix, denoted A, of (12) 
is full-rank, then (12) has a unique equilibrium. De¬ 
note this equilibrium xq = [wjf, Vq^, r/]]’, 0 q^]^. De¬ 
fine X = [uj"'- = [w^, v"^) ~ 

[u:l,VY,vlA'eV- Hence, 

X = Ax, (13) 

with the origin as its unique equilibrium. We are now ready 
to show the main stability result of this section. We first 
make the following assumptions. 

Assumption 1. The Laplacian matrix satisfies £0 = k^pEn. 



Assumption 1 can be interpreted as the emulated AC 
dynamics of (7) having the same susceptance ratios as the 
conductance ratios of the HVDC lines. 

Assumption 2. The gain 7 satisfies 7 > ■ 

Assumption 2 lower bounds for the damping coefficient 7 . 
Theorem 1. If Assumptions 1 and 2 hold, the origin of (13) 
is globally asymptotically stable. 


Proof. Consider the Lyapynov function candidate: 

1 p^nom 




(14) 


where C = diag(C'i,..., C„). Clearly W{uj,V is 
positive definite and radially unbounded. Differentiating 
(14) along trajectories of (13), we obtain 

w(cd, V, f], 0 ) = 

+ fj^fj + 

= uF{- + if'^™°P)a; 

+ - ifdroop. I - _ 

+ - {V^°^Lr+K'^)V + 

+ j^T(^^droop. 

+ $^S'^C4,s{s'^{K^)-^K‘^uj - 7 (^) 

= -uF{ - + A:‘i''°°P)w 

+ 2w^X“y - V'^{V’^°’^Cr + K^)V 
+ - fi'^Cn'h - 

since SS'^ = In - ^Inxn, so £ 4 , 33 '^ = C 4 ,. By defining 


1 ^', 


we obtain VJCrV =_ V"^3^Cr3V"_ and £ 4 , 34 ) = 
V"^3^£^3(j), where V" = [^ 2 ', • • •, By invoking 

Assumption 1 , we obtain 

w{£,v,nA) 



— 1 ^’ 



' 1 

V' = 

7®. 
_ 1 

V 

V = 

—^In 3 
_Cn 





CO 


y 




- v'^£r,V- 


=Qi 

V^°^3^£r3 -^3^£r3 
-^3^£r 3 jk^3^£R3 


V" 


= Q 2 


By applying the Schur complement condition for positive 
definiteness, we see that Qi is positive definite iff 


+ ^droop) _ k‘^[Kv)~^K‘^ 

= > 0 . 

Hence Qi is positive definite. By applying the same 
argument to Q 2 , Q 2 is positive definite iff 


Clearly the above matrix inequality holds under Assump¬ 
tion 2, since 3^£r3 > 0, and 3x ^ kin for fc 7 ^ 0. Thus 
W {ui, V, ry, ^) < 0, and the set where IT (a), V, fj, $) is non¬ 
decreasing is given by 

G = {(w, T,7,0)|lT(w,T,r7,^) = 0} 

= {(w, T,7,^)|i7 = 

for any fc € M. Clearly the largest invariant set in G is the 
origin, and thus fc = 0. By LaSalle’s theorem, the origin is 
globally asymptotically stable under the dynamics (13). 

We now turn our attention to the equilibrium of (12). 
If the stability condition in Theorem 1 is met, then the 
equilibrium is stable, and A is thus Hurwitz. Hence (12) 
must have a unique equilibrium. In the following corollary 
we show that the equilibrium of (12) satisfies Objective 1. 
Corollary 2. Let Assumption 1 hold and let 7 , k 4 , be 
given such that Assumption 2 holds. Let K^,K‘^ and 

^droop ^Yiat = iyp'(X‘^)-iiydroop 

, where F^ = diag{f[’,...,f^) and F^ = 
diag(/X ,... ,fX)- Then Objective 1 is satisfied in the limit 
when |[(iL‘^)“^ AT'^ 11 ^ —>■ 0 . 

Proof. By Theorem 1, (12) has a unique and stable 
equilibrium. The last n — 1 rows of the equilibrium imply 

= 0n-i. 

Now II(Ar‘^)“^AT'^ 11^ —> 0 in the above equation implies 
that 3'^io = 0 ^ (jj = kiln for some ki G M. Consider the 
(2n -I- l)th to 3nth rows of the equilibrium of (12): 

if droop, = 

By inserting to = kiln and premultiplying the above 
equation with 1 ^, we obtain that ki = 0 , so Cj = 0 „ and 
Equation (3) of Objective 1 is satisfied. This implies that 
V = k 2 ln for some G IR- Finally we consider the (n-l-l)th 
to 2 nth rows of the equilibrium of ( 12 ): 

- (£r + V + -^£ 4 , 34 )" = On (15) 

^nom y ^ ^nom y p^nom ‘r il \ / 

Inserting w = 0„ and premultiplying (15) with 1^ yield 
l^A:^f = 0 . (16) 

Inserting w = 0„ and 7 = feln in ( 6 ) yields 

pg®" = -A: 2 A:'^(A:“)-iA:'^™°P’il„, (17) 

where ps®" = [Pf®",..., P®®"]^. It now remains to show 
that the equilibrium of ( 12 ) minimizes the cost functions 
(4) and (5). Consider first (4), with the constraints Pf*^"-|- 
_ ^inj ^ 0 W = l,...,n and Er=i = 0 ;^ By 
summing the first constraints we obtain Ei=i = 

- Er=i -fT- The KKT condition of (4) is 

pPpgen ^ 

Since (F^)"^ = (17) and (18) are iden¬ 

tical for ^2 = kg. We conclude that (4) is minimized. Since 
pg®n = -K^{K‘^)-^K’^'^°°Pri = -A: 2 F'^(A:“)-iA:‘i™°Pl„ 
and w = 0 „, premultiplying the first n rows of the equi¬ 
librium of (12) with M“^, and adding to the (n -I- l)th to 
2 nth rows premultiplied with yields 

-v'^°'^£rV - fc 2 A:^(A:"’)”^A:'^™°pi„ = p“. 



Premultiplying the above equation with 1^ yields k 2 = 
- Er=i Pr E”=i j^v^^droop ■ Additionally, CrV is uniquely 

determined. Now consider (5). Note that and hence 
are uniquely determined by (4). By the equilibrium of 
(1), CrV = /“j, where P"j = ,..., . Thus, the 

KKT condition of (5) is 

F^V = Crt, (19) 


in the simulation also consider the inductances and 
capacitances Cij of the HVDC lines. The capacitances of 
the terminals are assumed to be given by Ci = 0.375 x 10“^ 
p.u., and are chosen uniformly for all VSC stations. The 
AC grid parameters were obtained from [Milano 2010]. The 
generators are modeled as a 6th order machine model con¬ 
trolled by an automatic voltage controller and a governor 
[Kundur 1994]. Each load in the grid is assumed to be 
equipped with an ideal power controller. 


where r € K". Since CrV is uniquely determined, we pre¬ 
multiply (19) with 1^ and obtain the equivalent condition 

l^F^V = 0. (20) 

Since F^ = (16) and (20) are equivalent. Hence (5) 

is minimized, so Objective 1 is satisfied. 

Remark 4- Corollary 2 provides insight in choosing the 
controller gains of (6) and (7), to satisfy Objective 1. 

4. SIMULATIONS 

In this section, simulations are conducted on a test system 
to validate the performance of the proposed controllers. 
The simulation was performed in Matlab, using a dy¬ 
namic phasor approach based on [Demiray 2008]. The 
test system is illustrated in Figure 1. The line parameters 



Figure I. MTDC test system, consisting of a 6-terminal 
MTDC grid. Each terminal is connected to an IEEE 
14 bus AC grid, sketched as octagons. 

of the MTDC grid are given in Table 1. Note that we 


Table 1. HVDC grid line parameters 


i j Rjj [p-n.] Lj 
1 2 0.0586 

1 3 0.0586 

2 3 0.0878 

2 4 0.0586 

2 5 0.0732 

2 6 0.1464 

3 4 0.0586 

3 5 0.1464 

4 5 0.0732 

5 6 0.1464 


[10 3 p.u,] Cjj [p.u.] 


0.2560 

0.0085 

0.2560 

0.0085 

0.3840 

0.0127 

0.2560 

0.0085 

0.3200 

0.0106 

0.6400 

0.0212 

0.2560 

0.0085 

0.6400 

0.0212 

0.3200 

0.0106 

0.6400 

0.0212 


Table 2. Controller parameters 


T^LJ T^V ^droop 

1000 Too 9 


^droop, I 
_ 

3.35 


c”. 


5/Ri, 


4> 

C-- 
L2 


15/Pi. 


7 

0 


The test grid was controlled with the controllers (6) and 
(7), with parameters given in Table 2. The communication 
network of (6) and (7) is illustrated by the dashed lines in 
Figure 1. Note that we have set 7 = 0, and that Theorem 1 
thus does not guarantee stability of the equilibrium. The 
resulting matrix of the closed-loop system is however veri¬ 
fied to be Hurwitz. At time t = 1 the output of one gener¬ 
ator in area 1 was reduced by 0.2 p.u., simulating a fault. 
Figure 2 shows the average frequencies of the AC grids, 
while Figure 3 shows the DC voltages of the terminals. 
Figure 4 shows the total increase of the generated power 
within each AC area. It can be noted that immediately 
after the fault, the average frequency of the AC area of the 
fault drops. The frequency drop is followed by a voltage 
drop in all terminals, and a subsequent frequency drop in 
the remaining AC areas. The frequencies converge to the 
nominal frequency, while the voltages converge to their 
new stationary values after approximately 30 s. We note 
that the frequencies are restored to the nominal frequency, 
as predicted by Corollary 2. Furthermore the incremental 
generated power of the AC areas converge to the same 
value, as a consequence of Corollary 2 and the controller 
parameters being equal. 


5. DISCUSSION AND CONCLUSIONS 


In this paper we have studied a distributed secondary 
controller for sharing frequency control reserves of asyn¬ 
chronous AC systems connected through an MTDC sys¬ 
tem. Under certain conditions, the proposed controller 
stabilizes the closed-loop system consisting of the intercon¬ 
nected AC systems and the MTDC grid. The frequencies 
in all AC grids are shown to converge to the nominal 
frequency. Furthermore, quadratic cost functions of the 
voltage deviations of the MTDC terminals and of the 
generated power, are minimized asymptotically. Finally, 
the results were validated on a six-terminal MTDC system 
with connected AC systems. 











































Figure 2. Average frequencies in the AC areas. 



Figure 3. DC terminal voltages. 



Figure 4. Increase of generated power in the AC areas. 
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